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Abstract 



> 

m 

By a Wigner-function calculation, we evaluate the trace of a certain Gaussian op- 
erator arising in the theory of a boson system subject to both finite temperature 
£Sj . and (weak) interaction. Thereby we rederive (and generalize) a recent result by 

Kocharovsky, Kocharovsky, and Scully [Phys. Rev. A 61, 053606 (2000)] in a way 
that is technically much simpler. One step uses a special case of the response of 
Wigner functions to linear transformations, and we demonstrate the general case 
by simple means. As an application we extract the counting statistics for each mode 

' of the Bose gas. 
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In a recent paper on dilute boson gases [1], Kocharovsky, Kocharovsky, 
and Scully faced the problem of evaluating the trace of a product of two 2- 
mode Gaussian operators. They solved it by exploiting the algebra of 2-mode 
squeezing operators in a remarkable tour de force occupying the appendix of 
[1]. The main objective of the present Brief Report is to put on record an 
alternative derivation that is simpler and more direct. 

In a few words, the essence of our calculation is this: We recall (i) that the 
trace of a product of two operators can be obtained as a phase space integral 
of the product of the respective Wigner functions [2]; (ii) that the Wigner 
functions of Gaussian operator functions are Gaussian functions of the phase 
space variables; and (iii) that integrals of Gaussians are immediate. 
The trace of interest is 

T(u, v- e) = tr{GK v)p e } = (G(u, v)) e (1) 
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Table 1 

Dictionary for switching from the present notation to the one in Eqs. (63), (68), 
(69) of Ref. [1]. 



here 


there 


here 


there 


T(u, u; e) 


e± k («) 


u, V 


u, u 


a, at 
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£k/T 


6, &t 


/3-k,/?l k 


tanh$ 


A k 


A, A* 




Y+ 


z(-A k ) 


B, 


&-k,&l k 


Y_ 





with 

G(u,v) =e iuat ° + iv6t6 (2) 

and 

Pe = (l-e-*) 2 e-< A]A + B * B ), (3) 

where at , a and 6 are the ladder operators of the physical excitations (two 
independent harmonic oscillators), and the Bogolubov transformation 

A = a cosh •& — sinh # , A* = at cosh •& — b sinh $ , 

B = b cosh ■& - at sinh # , S 1 " = 6 1 " cosh ft - a sinh $ (4) 

(with real d) relates them to the ladder operators of the "dressed modes." The 
statistical operator p e is a thermal state (with temperature parameter e > 0) 
that is diagonal in the dressed number operators A* A and EtB. Table 1 states 
how our notation is related to the one used in [1]. The factor G(u,v) is the 
generating function for the count of physical excitations inasmuch as 

^a,Jm, n - j ^ j ^-' m ~' mV GM, (5) 

(2tt) (2tt) 

the integrations covering any interval of length 2n. Accordingly, upon expand- 
ing T{u, v; e) in powers of e m and e lv we get the probabilities for having certain 
numbers of quanta in the physical excitations. 

1 Phase-space evaluation 

Our starting point is the 1-mode thermal state, 

p?) = (1 - e- e )e-- ta , (6) 
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and its well known Wigner function, 

[p«] w = 2tanh(ie)e- 2a * atanh ^ e ) , 
normalized in accordance with 

d 2 a 



tr{p«}=/^%, 

J 71 



(1)1 

€ JW ' 



(7) 



(8) 



where a is the complex phase space variable associated with a, and the stan- 
dard parameterization of the phase space integral — specified by a = x + iy, 
a* = x — iy, and d 2 a = dx dy — suffices for our purposes. 

First, we need the Wigner function (or Weyl transform) of G(u,v) in (2), 
which (among other methods) can be found efficiently by analytic continuation 
from (6) and (7) [taking e — > — iu and e — > — iv in two copies of (7), discarding 
its prefactor]: 



[G(u, v)] w = e-^ u + v ) sec(^) sec(^) 

xe 2ia*atan(iu) + 2i(3* (3 t&n(±v) 



(9) 



where (3 is the phase space variable for b. 

Second, to get the Wigner function of the statistical operator (3) we recall 
a fundamental property of Wigner functions, namely that linear similarity 
transformations of the ladder operators carry over to the Wigner function 
(see Sec. 4 below). For the Bogolubov transformation (4) this means that 
[F(at,o,6t,6)] w = f( a *,a,P*,P) implies [3] 



[F{A\ A, B\ B)] w = f(a* cosh 0-0 sinh #,...,/? cosh - a* sinh ■&) , 



(10) 



and therefore we have 



[p e ] w = 4tanh 2 (|e) 



xe 



-2[(a*a + (3* (3) cosh(2tf) - (a*f3* + af3) sinh(2t?)] tanh(±e) 



Given (9) and (11), the trace (1) is a 4-dimensional Gaussian integral that 
is routinely evaluated, 



T(u,v;e) 



d 2 a r d 2 /3 



n 



n 



[G(u,v)] w [p t ] w 



= 4(e f 



(e e -l) 2 (l + e m )(l + 
+ (e e + 1) 2 (1 -e iu )(l-e [v ) 



+ 2(e 2e - l)(l-e 



iu + iv 



) cosh(2tf) 



-i 



(12) 
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While this formula is the most direct and natural result of the calculation, it 
may appear physically inscrutable. However, it can be worked into the form 

rp/ v 2(y + -i)(y_-i) 

T(u,v;e) ~ 



(Y+ - e m ){Y^ - e w ) + (y_ - e m )(Y + - e w ) 
= T{v,u;e), (13) 

where 

e e + tanh-i9 e e — tanht? 

y — y — (14) 

+ l + e e tanhtf' " 1 - e e tanhtf ' V ; 

This matches the result, for the special case u — v, reported in (69) of [1] 
(translated with Table 1): 

The challenge of rederiving (15) by simpler means has thus been met; in fact, 
(13) is more general because u and v need not be the same. 

Remark: One can also calculate the trace (1) by using the Q-function for 
p e and the P-function for G(u,v), rather than the Wigner function for both. 
But that approach is more involved because the statement corresponding to 
(10) is more complicated. 



2 Occupation probabilities 

Let us now see what we can learn from this gain in flexibility. Special 
cases are those of v = and u — 0, the generating functions for the count of 
excitations in either the a mode or the b mode, 

%0; £ ) = (e Wfl ) e V 1 
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T(0,v,e) = (^ b \ = ^— I, (16) 



where 



Y - 2Y+Y - ~ Y +- Y - e e + tanh 2 g 

Y + + Y_-2 l+e e tanh 2 ^ ' 1 } 

The symmetry in (16) is, of course, just a particular example of T(u,v;e) = 
T(v,u; e), which we noted at (13) in passing. The geometric series that result 
from expanding the right-hand sides of (16) in powers of e m or e lv , respectively, 
tell us the individual counting statistics for each mode, 



^ta,n) £ = ^t 6 , n ) £ =(l-^)^. (18) 
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So, if there were no statistical correlations between the two modes, the prob- 
ability for having a total of N quanta in both modes would be given by 



N 

Wja, „>«%(,, *-„>. = <*+ DC - Y- 1 )^-" ' < 19 ) 

n=0 



But there are correlations and, indeed, the correct value of {^ a \ a _|_ ^5 
obtained by expanding T(u,u; e) of (15) in powers of e m , 

<*ata + 6t6, iV>, = ^ " ^ " ^) ^^^ , (20) 

differs noticeably from (19), unless we are in the limiting situation of 1? = 
when Y + = Y_ = Y. In that limit we have 

T(u, v- e)\# = = f ~ a f ~ I , (21) 

and this factorization states the obvious: For $ = the modes do not get 
dressed and so they maintain their statistical independence. 

Rather more involved is the expression for the probability of m quanta in 
mode a and n quanta in mode b, 



~ (m + n-k)\ 

x £ 



^ fc!(m - k)\(n - k)\ [(Y + + Y_f\ ' 



l^m+n 
I k 



(22) 



where the summation terminates when k equals the smaller one of m and 
n. Note that the right-hand side is invariant under the interchange m <-> n, 
which is another manifestation of T(w, v; e) = T(i>, it; e). Incidentally, one could 
write (22) in terms of a Jacobi polynomial, but that doesn't seem to add 
transparency to the result. 



3 Bogolubov transformation of a single mode 

Also worthy of study is a Bogolubov transformation mixing one mode with 
itself by 

A = a cosh -d — a) sinh d (23) 

in place of (4). In fact, the dressing of a pair of modes as in (4) can be di- 
agonalized into two decoupled transformations of the form (23); in the Bose 
condensation problem this corresponds to replacing traveling waves by stand- 
ing waves. 
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In this case we should evaluate 

T(u;e)= [^[G(u)] w [p e ] w , (24) 

J 7T 

where is the v — version of (9) and 

[ Pt ] w = 2tanh(|6)e-[ 2 «*« cosh (2^) " («* 2 + « 2 ) sinh(2^)] tanh(^) (25) 

is obtained by a — > a cosh?? — o;*sinh$ in (7). [Note the subtle difference 
in the relative sizes of the terms in the exponent in (25) versus (11).] The 
Gaussian integration is more elementary than in the previous case and yields 

T(u; e) = 2(e e - 1) [(e e - 1) 2 (1 + e iu ) 2 + (e e + 1) 2 (1 - e m f 



+ 2(e 2e -l)(l-e 2iw )cosh(2tf)j 2 (26) 

as the analog of (12). Note that T(u; e) 2 = T(u, u; e), as it should. On the other 
hand, T(u; e)T(v; e) does not equal T(u, v; e), as these functions are giving the 
statistics for two different bases of single-quantum states, which have different 
dynamics. If desired, the counterpart of (22) could be obtained by expanding 
(26) in powers of e m , which is easily done because we meet here a familiar 
generating function for Legendre's polynomials. 

One motivation for investigating a single mode was the hope of elucidating 
the structure of equations (14), which cry out to be interpreted as addition 
formulas for the hyperbolic tangent function with e e = tanh(. . .). There is, 
in fact, a way of constructing the thermal state of a boson mode by means 
of another Bogolubov transformation [4]; one would think that (14) then ex- 
presses the composition of the two transformations. However, the appropriate 
thermal parameter in that construction is \e rather than e. This discrepancy 
may be related to the fact, stressed in [1], that quanta in the interacting Bose 
gas are present in strongly correlated pairs. Calculations following [4] are in 
progress, but so far have not yielded expressions as simple as those provided 
by the method reported here. 



4 Linear similarity transformations 

It is striking that the fundamental transformation property of Wigner func- 
tions, of which (10) is the special case we make use of, is not mentioned in any 
of the reviews [5-7] or textbooks [8,9], although it has been much exploited 
for various applications (e.g., in [10,11]). Since it deserves to be known more 
widely, we present the following compact derivation, which is more general 
than all previous ones. 

As a warm-up we'll first consider the single- mode case of Sec. 3. The Wigner 
function (or Weyl transform) F w (a*,a) of an operator function F(a\a) is 
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given by 

F w (c**, a) = J e^* ~ F°ti {eP a " & F{a\ a)) 

= 2tr|e aat - a * a 5e a * a - aat F(at,a)| (27) 

* _ t 

where we encounter the unitary displacement operator e a a aa , 

Q aa) — a* a a \ e a*a — aa^ _ Q f _ a * 

e aat - a *a ae a*a-aa^ = a _ a ^ ( 2 8) 
and the reflection operator 

S= (*l<Jr«-P«\ (29) 
J 2n 

which is both unitary and Hermitean and therefore its own inverse, & = S = 
S^ 1 [12]. Its name derives from the unitary transformation that S effects, 

S- 1 atS = -at, S^aS^-a. (30) 

A linear similarity transformation turns F(a),a) into G(a\a), another 
operator function, 

G(a\ a) = V~ 1 F{a\ a)V = F(fj,a) + va, aa + ra f ) , (31) 

where 

V~ x aW = fia 1 + va , VafV -1 = aa) — fa , 
V' l aV = aa + ra^ , VaV' 1 = fia - ra^ . (32) 

The commutation relation [a, a)] = 1 puts the restriction /ja — vt = 1 on the 
numerical coefficients fi, v, a, r, and otherwise they can take on any complex 
values permitted by \a\ > \v\, |r|, which condition ensures that V~ x a)V has 
eigenbras and V~ 1 aV has eigenkets. It is not necessary that (32) be a unitary 
transformation, such as (23). In particular, the non-unitary transformation 
a) — > a) — a, a — > |(a + a)) is permissible and useful [10]. 

The crucial observation is now that the reflection operator is invariant 
under this similarity transformation, 

V -±SV = S or S- 1 VS = V. (33) 

This is so because V is the exponential of a linear combination of a) 2 ', a)a, 
and a 2 (or a product of such exponentials), as illustrated by 

V = e ^ r /^ at V~ a1a e~^W^ a2 . (34) 
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Clearly, the reflection (30) has no effect on V, which is the second statement 
of (33). 

Now, perform the replacements 



F(a\ a) -> V~ 1 F{a\ a)V = G(at, a) , 

s -> s = v- l sv 



(35) 



in (27) to establish 

G w (a*,a) = 2tr^Ve aat - a * a V- 1 SVe a * a - aa W- 1 F(a\a)^ , (36) 



then note that 



y e a * a ~ aa^ y-i _ e (/ua* + va)a — [pa + ra*)a) 



and arrive at 



Gw(«*, a) = Fw(fj,a* + ua, aa + to*) 



(37) 



(38) 



Indeed, the linear similarity transformation (31) of the operator functions 
carries over to their Wigner functions. 

The arguments in higher dimensions are quite analogous. There are now n 
pairs of ladder operators and n pairs of phase-space variables associated with 
them. We write them compactly as 2n-component rows, 



a = {a\,ai,a\,a 2 ,...,a\ l ,a n ) 



a = (al, a u a* 2 , a 2 , . . . , a*, a n ) . 
The n-dimensional displacement and reflection operators are then 

7J e a l a k ~ a k a\ = g aKa T ^ S= f 

k=i •* (2n) n 

with the 2n x In block matrix 

'0-1 
1 



(39) 



(40) 



K 



V 



-1 

1 o 



-K- 1 = -K 1 



(41) 



so that the analog of (27) reads 

F w (a) = 2 n tr | e - 5K « T 5e 5K « T F(a)| 



(42) 
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A linear similarity transformation in n dimensions is of the form 

V- x aV = a V , VaV' 1 = a V" 1 (43) 
where V is a 2n x 2n matrix restricted mainly by 

VKV T = K , (44) 

which is the analog of /icr — vr = 1 for the coefficients in (32). Here, too, V 
is the exponential of a bilinear form of the ladder operators, V = exp(aVa T ), 
with a (largely) arbitrary complex 2n x 2n coefficient matrix V, and so (33) 
is equally true in n dimensions. By arguments analogous to those of (35)-(38) 
it then follows that the Wigner function of the transformed operator 

G(a) = V^F^V = F{a V) (45) 

is given by 

G w (a) = 2 n tr{ve- aKsT V- 1 SVe SKsT V- 1 F(a)} 
= 2™tr je-« VK « T Se 5VK « T F(a)j 

= F W («V) . (46) 

This closes the case. 
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